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From Planar Solitons to Vortex Rings and Lines: 
Cascade of Solitonic Excitations in a Superfluid Fermi Gas 
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We follow the time evolution of a superfluid Fermi gas of resonantly interacting ®Li atoms after 
a phase imprint. Via tomographic imaging, we observe the formation of a planar dark soliton, its 
subsequent snaking, and its decay into a vortex ring, which in turn breaks to Anally leave behind 
a single solitonic vortex. In intermediate stages we And evidence for an exotic structure resembling 
the 4>-soliton, a combination of a vortex ring and a vortex line. Direct imaging of the nodal surface 
reveals its undulation dynamics and its decay via the puncture of the initial soliton plane. The 
observed evolution of the nodal surface represents dynamics beyond superfluid hydrodynamics, 
calling for a microscopic description of unitary fermionic superfluids out of equilibrium. 


Solitonic excitations such as solitons, vortices and vor¬ 
tex rings are found in a large variety of nonlinear media, 
from classical fluids and plasmas to polyacetylene chains 
and superconductors. While ubiquitous, their intrinsic 
properties are tailored by the host medium. In superflu¬ 
ids, which are characterized by a complex order parame¬ 
ter with a well-defined phase and a non-viscous flow, such 
excitations correspond to phase defects and exhibit prop¬ 
erties non-existent in their classical counterparts. There, 
a vortex is topologically protected owing to the quantized 
circulation of the velocity field, and a traveling soliton ex¬ 
periences superfluid back flow determined by the phase 
difference across it mm. The quantum statistics of the 
particles forming the superfluid is yet another ingredient 
which dramatically affects the properties of these defects. 
In Fermi superfluids, as opposed to the bosonic case, dark 
solitons and vortices are known to host in-gap fermionic 
excitations in their cores, from the Andreev bound states 
in the generic case mm, to the more exotic Major ana 
fermions in the presence of spin-orbit coupling EE]. 

Importantly, in a quantum fluid with short-ranged in¬ 
teractions, these phase defects are localized within the 
microscopic length scale of the system: the healing length 
f. The healing length sets the length scale above which 
the superfluid dynamics is well captured by the hydro- 
dynamic formalism. At length scales on the order of f or 
smaller, a microscopic description is required, and this is 
where the dichotomy between Bose superfluids and Fermi 
superfluids becomes stringent. While weakly interacting 
Bose-Einstein condensates (BEC) are well understood in 
terms of the Gross-Pitaevskii (GP) theory, a complete 
microscopic wave equation for strongly-interacting Fermi 
superfluids remains to be established. At the mean- 
field level, a unified description can be formulated within 
the Bogoliubov-de Gennes (BdG) formalism, which con¬ 
nects to the GP equation in the limit of weakly inter¬ 
acting BEGs, and contains the necessary fermionic de¬ 
grees of freedom in the Bardeen-Gooper-Schrieffer limit 
(BGS) dJ E E [7]. However, while the BdG framework 
provides a good description of these two limiting cases, 


it is unclear whether it contains the right ingredients to 
quantitatively handle the behavior of solitonic excitations 
in the strongly correlated regime, where the dynamics 
near the core of these phase defects is highly nontriv¬ 
ial [8]. The unitary Eermi gas realized in ultracold atom 
experiments offers a unique opportunity to clarify this is¬ 
sue, as it resides at the point of the BEG-BGS crossover 
where beyond mean-field correlations are expected to be 
the strongest [9]. It is also the regime where the healing 
length f is the smallest - on the order of the inter-particle 
spacing - such that phase defects are as localized as pos¬ 
sible in a quantum fluid. 

A natural approach to experimentally reveal the core 
dynamics of such defects is to trigger their decay. Soli¬ 
tonic excitations indeed follow a well-defined hierarchy 
in terms of stability and energy cost in three dimensions, 
the planar soliton being the most energetic and unstable 
towards the formation of other solitary waves pHUTb] . 
In weakly interacting BEGs, dark solitons have been ob¬ 
served to decay into vortex rings and vortices pIMTS] as 
a consequence of the snake instability, the undulation of 
the soliton plane m- In the case of strongly interacting 
Fermi superfluids, similar scenarios have been predicted 
numerically within a mean-field approximation [MU], 
but an experimental support of such microscopic dynam¬ 
ics is still lacking. 

In this Letter, we create a cascade of solitonic excita¬ 
tions in a unitary Fermi gas of ^Li atoms. Starting from 
a planar dark soliton created via phase imprinting, we 
observe the formation of ring defects which eventually 
decay into a single solitonic vortex. By means of a to¬ 
mographic imaging technique [22], we are able to follow 
the surface dynamics of the soliton’s nodal plane at the 
level of the interparticle spacing, as it snakes, breaks and 
converts into the topologically protected solitonic vortex. 
Our measurements allow for a quantitative analysis of 
the snaking dynamics of the initial dark soliton, awaiting 
comparison to time-dependent theories of strongly corre¬ 
lated fermions. 

We create a strongly interacting fermionic superfluid 
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FIG. 1. (Color online) (a) Cascade of solitonic excitations in a unitary Fermi superfluid following the phase imprint. A planar 
soliton snakes and decays into a vortex ring. Shown are images of the density distribution in the central slice of the superfluid, 
after rapid ramp and time of flight, for the first 20 ms after the imprint. The imprint also generates two sound waves propagating 
towards the edges, (b) Time-series of the central slice up to t = 100 ms, cropped to the region around z = 0. 


using a balanced mixture of the two lowest hyperfine 
states of ^Li (|1) and |2)) at a Feshbach resonance [23] . 
Our atomic clouds contain 7 x 10^ atoms per spin state 
confined in an elongated trap, combining a tight radial 
optical potential (in the x-y plane) and a shallower ax¬ 
ial magnetic potential (along the ^ axis). The axial and 
radial trapping frequencies are = 10.87(1) Hz and 

= 69(6) Hz, respectively. The axial and radial 
Thomas-Fermi radii of the cloud are Rz = 326(2) ym 
and R± = 54(2) /xm, and correspond to a chemical po¬ 
tential at the center of the cloud of y — hx 3.7(1) kHz = 
54(4) hLJ±. The gas is thus deep in the three-dimensional 
regime. Gravity slightly weakens the trapping potential 
along the vertical ^-direction (see [24]). Phase imprinting 
is realized as in Refs. [2^1251128] . whereby one half of the 
superfluid is exposed to a blue-detuned laser beam for 
a time sufficient to advance the phase by approximately 
TT. These experimental parameters are similar to those 
of previous works [22 EH], where a single solitonic vor¬ 
tex was detected and observed to undergo a determin¬ 
istic precessional motion for several seconds. Here, we 
study the evolution of the excitations at early times fol¬ 
lowing the phase imprint. To probe such dynamics, we 
employ the detection scheme introduced in [22], which 
combines the so-called rapid ramp technique and tomo¬ 
graphic imaging [24]. The rapid ramp converts the spa¬ 
tial variations of the order parameter A(r) of the fermion 
pair superfiuid into resolvable density variations of a Bose 
gas of molecules [28fi3Q] . It is a powerful tool for the 
study of phase defects in atomic Fermi gases, as the 
presence of any phase singularity manifests itself as a 
strong density depletion. In our experimental sequence, 
the rapid ramp is performed at a variable wait time t 
following the imprint. We then slice a thin layer of the 


atomic cloud at a chosen y position, and destructively 
probe its density distribution via absorption imaging. 

Fig.[TJi shows a time sequence of images recorded in the 
first 20 ms after the phase imprint, which correspond to 
the density distribution at the central slice (near y — 0) 
of the superfiuid. At the location of the imprinted phase 
jump ( 2 ^ = 0), a slow and straight dark soliton emerges 
and subsequently undergoes a snaking motion, seeding 
the puncture of the nodal surface. The broken soliton 
evolves into a vortex ring structure, visible in the cen¬ 
tral slice as a pair of nodal points. Fig. presents 
a zoomed-in view of the soliton’s time evolution up to 
100 ms after the imprint. Simultaneous with the soli¬ 
ton’s core dynamics, two wavefronts quickly propagate 
to the edges of the cloud, which we identify as sound 
waves. The upper and lower sound wavefronts are found 
to propagate at speeds of 13.1(4) mm/s and 13.1(8) mm/s 
respectively, which coincide with the speed of sound of 
12.9(1) mm/s estimated from the peak density using the 
relation Cg = where .^b = 0.37 is the Bertsch 

parameter [31] and v-p the Fermi velocity [32] . The appar¬ 
ent large amplitude of these sound waves is a consequence 
of the rapid-ramp M- 

The dynamics is analyzed in detail in Fig. [^ show¬ 
ing residuals of the central slice as a function of time, 
along the axial cut at x = 0 (Fig. If) and along its 
outer edge near x — R± (Fig. §) . The different char¬ 
acteristic speeds of the various waves generated after the 
phase imprint are apparent. One recognizes the two ini¬ 
tially created sound waves following linear trajectories 
with opposite slopes, while the dark soliton remains near 
z = 0 with negligible velocity. A second set of shallower 
sound waves is emitted about ~ 5 ms after the initial 
sound wavefronts, forming all together a pattern of hy- 











3 



FIG. 2. (Color online) Overview of the dynamics following 
the phase imprinting, (a) Representative images of the central 
slice at t = 4, 50 and 95 ms showing a planar soliton (S), the 
nodal points of a vortex ring (R) and a vortex line (V). The 
vertical lines indicate the regions of interest of the generation 
of residuals in (b). (b) Residuals on the central slice along the 
central axial cut x = 0 (upper panel) and an outer axial cut 
(lower panel). Dark (bright) color indicates density depletion 
(excess). Two sound waves rapidly propagate to the edges, 
while a sharp depletion remains at the center. Around t = 
5 ms, a second set of shallower sound waves is emitted. The 
residuals show the puncture of the soliton plane in the central 
slice (t ~ 15 ms) and the return of a vortex line at t ~ 80 ms. 


drodynamic wakes. The rapid vanishing of the sound 
contrasts with the persistence of the solitonic wave near 
2 ; = 0. The unbroken soliton appears in the residual near 
X = 0 (Fig. Hf) as well as the outer (Fig. H) residuals, 
while the vortex ring structure is signaled by a localized 
point-like depletion visible exclusively in the outer resid¬ 
uals. Close to 100 ms after the imprint, a single solitonic 
vortex remains, precessing in the superfluid, with a pe¬ 
riod of ~ 1.4 s along the 2 ^ axis [22l [28] (see also EH), 
which is more than an order of magnitude longer than 
the duration of the cascade. 

In order to obtain a complete picture of the dynamics, 
three-dimensional tomography of the superfluid is per¬ 
formed after the rapid ramp. Fig. displays a set of 
representative tomographic images at various times, giv¬ 
ing access to the local pair density after the ramp. From 
these images, we are able to reconstruct the structure of 
the defect engraved in the superfluid and follow its time 
evolution. The right panel of Fig. shows the recon¬ 
structed depletion as would be seen from the long axis of 
the cloud. At early times, a surface of depletion cutting 
through the entire cloud’s section is observed, the planar 
dark soliton. It subsequently tears in its upper half and 
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FIG. 3. (Golor online) Tomography of the cascade. Main 
panel (left): representative tomographic images at different 
stages of the cascade. Right column: structure of the de¬ 
pletion due to the defect, as would be seen along the z-axis, 
reconstructed from the tomographic images, t = 3 ms: sharp 
density depletion across the whole cloud signaling a planar 
dark soliton. t = 7 ms: snaking of the soliton plane and first 
signature of a puncture, t = 9, 16 and 20 ms: the punc¬ 
ture in the upper half of the soliton plane grows and yields 
an asymmetric vortex ring, with a nodal area left at bottom, 
t = 85 ms: the lower nodal area is punctured and a vortex line 
forms across a vortex ring, t = 95 ms and t = 400 ms: the ring 
part of the defect progressively disintegrates, t = 1000 ms: a 
single solitonic vortex remains. 



then undergoes a cascade into structures with smaller 
and smaller nodal area. The hole appearing in the initial 
nodal surface is seen to continuously grow in size, leading 
to the formation of a transient asymmetric vortex ring, 
combining the bottom part of the initial soliton and a 
vortex line bent in a semicircle on the upper part. At 
this stage one might anticipate that the nodal area will 
naturally heal into a standard vortex ring, by shrinking 
into a single loop with a core of size f. However, the to¬ 
mographic images obtained at later times suggest a more 
complicated scenario [24], where a second puncture oc¬ 
curs in the lower nodal plane. This results in the forma¬ 
tion of a horizontal line depletion which we interpret as 
a vortex line intersecting the vortex ring. This structure 
is seen in the tomographic images e.g. at t = 85 ms (see 
Fig. [3|and resembles the ^-soliton recently proposed in 
Ref. [15] • At even later times, the ring part of this exotic 
defect progressively disintegrates (as seen at t = 95 ms 
and 400 ms in Fig. [^, leaving behind a single solitonic 
vortex (t = Is), which precesses in the superfluid. It is 
the precession of this remnant solitonic vortex which has 
been studied in [221128]. 

Recently, several theoretical works have treated the 
evolution of fermionic superfluids following a phase im¬ 
print [m EH ESI ED and the possible cascade scenarios 
following the decay of a planar dark soliton [201I2T1I351I36] 
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via various mean-field approaches. In some of these 
works, it has been found numerically that in a cylin- 
drically symmetric potential with negligible dissipation, 
a planar soliton decays into a vortex ring, which then 
undergoes a long-lived oscillatory motion along the ^ 
axis [2S1 ISni EH]- By mimicking experimental imperfec¬ 
tions, such as trap distortion due to gravity uni, and 
imperfect phase imprinting [211 EH, later works found 
that the vortex ring further decays into a single rem¬ 
nant vortex. The proposed scenarios are however distinct 
from our observations. Recent simulations based on the 
Gross-Pitaevskii equation reveal a multitude of dynami¬ 
cal pathways towards the final single vortex, via various 
intermediate “Chladni solitons” [laEi]. 

At the origin of the cascade lies the snaking instabil¬ 
ity of the soliton. In order to quantify this undulation 
dynamics, we perform a Fourier analysis on the shape 
of the soliton Zs{x) ^ Aq cos 

in terms of the transverse modes of wavelengths = 
2R^/n, the integer n being the mode number, with 
Fourier amplitudes and phases 0^. Fig. displays 
selected images of the snaking soliton, and Fig. the 
corresponding nodal profiles ^s(^) obtained from the den¬ 
sity minima. For each of these curves, the result of the 
Fourier expansion up to the 5^^ order (gold solid line) is 
superimposed, illustrating that the undulation observed 
here is well characterized in terms of transverse mode 
excitations. Fourier spectra of the soliton’s undulation 
are obtained for 2 ms < t < Ilms (see Fig. |^), and the 
evolution of the amplitudes A^ is reported in Fig. Hi- 
We find that the fundamental mode Ai = 2R^ largely 
dominates this dynamics, with a relative weight consis¬ 
tently higher than that of the harmonics and a signifi¬ 
cantly larger growth rate Ai. The rates decrease as 
the mode number n increases (see Fig. Ef) . The contri¬ 
bution of the modes n > 5 was found to be insignificant. 

The trend of the growth rates (Fig. |^,e) is in 
strong contrast to what is expected for a uniform trans¬ 
verse confinement. In that case, a hydrodynamic ar¬ 
gument [38] implies that the growth rate at long wave¬ 
length increases linearly with the transverse mode [39] . 
as confirmed, for small wavevectors, in BdG numeri¬ 
cal simulations for crossover superfluids m- Physi¬ 
cally, the predominance of the fundamental mode found 
here is not surprising given the transverse inhomogene¬ 
ity of the trapped superfluid. Indeed, the phase differ¬ 
ence A(j) across the soliton is tied to its relative velocity 
Ug/c, where c is the superfluid critical velocity, through a 
unique current-phase relation A(j) = /(ug/c). While the 
functional form of / is still under debate for the unitary 
Fermi gas, it is expected to be monotonic, at least at 
small relative velocities Therefore, for an inho¬ 

mogeneous transverse confinement where c is highest at 
the center of the cloud (r = 0) and gradually decreases 
towards the edges (r ^ R±), the phase jump A(p across 




FIG. 4. (Color online) Spectral analysis of the snaking dy¬ 
namics. (a) Snapshots of the soliton’s undulation, (b) Ex¬ 
tracted undulation shapes Zs{x) (black dots) along with the 
corresponding Fourier expansions to the 5^^ order (solid gold 
line), (c) Fourier spectra at t = 3 and 11ms. (d) Fourier 
amplitudes An as a function of time, for n = 1 (circle), 2 
(square), 3 (triangle) and 4 (diamond) and their fits to a line 
(solid lines). The error bars indicate the standard deviation 
of the mean obtained from a 3-points binning of the data, (e) 
Growth rates An for n = 1 to 5, obtained from the linear fits 
in (d), with the error bars being the fit error. 


a moving soliton can only be maintained if its velocity 
is locally adjusted such that Ug(r)/c(r) is constant. The 
soliton thus develops a drum-like profile that continu¬ 
ously stretches along the propagation direction, which 
explains at least in part the growth of Ai. 

In conclusion, we have observed a cascade of solitonic 
excitations in a strongly-interacting Fermi superfluid, 
from an initial planar dark soliton towards a final, rem¬ 
nant solitonic vortex, through an intermediate ring struc¬ 
ture resembling the recently predicted 4>-soliton [15]. At 
the origin of the cascade lies the snaking instability, 
which we characterized quantitatively by studying the 
evolution of transverse Fourier modes. The breaking dy¬ 
namics of the unitary Fermi gas studied here occurs at 
the scale of the interparticle spacing, and provides a novel 
experimental input for microscopic theories of strongly 
interacting fermions. A natural extension of our work 
is to approach a regime where the snake instability is 
inhibited, e.g. via a strong confining potential in the 
radial direction. Future prospects are a measurement 
of the soliton’s current-phase relation in the BEG-BGS 
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crossover m, the detection and manipulation of An¬ 
dreev bound states trapped inside the soliton [ 2111 ], and 
the creation of soliton trains in the presence of spin im¬ 
balance, which would realize one limit of Fulde-Ferrell- 
Larkin-Ovchinnikov states [42H4^ . 
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SUPPLEMENTAL MATERIAL 


Supplementary data. In Figs. i 0 and we 

present supplementary data that could not be included 
in the main text. Fig. shows the time evolution of 
the central slice density, after rapid ramp and time 
of flight, from t = 0 to t = 19.5 ms in time steps of 
0.5 ms. In Fig. we show all the images used in the 
spectral analysis of the snake instability, along with 
their detected minima and Fourier reconstruction. In 
Fig-i we show the full tomography recorded at several 
instances of the dynamics from t — 3.5 ms to t = 1000 ms. 

^-soliton or not? The tomographic images taken at 
the instants t ~ 50 — 100 ms show both the presence of 
horizontal vortex lines and pairs of nodal points. Be¬ 
cause of the destructive nature of our measurement, one 
could a priori attribute this observation to two distinct 
situations: (A) the imprinting leads to either a vortex 
ring or a vortex line with unknown probabilities; (B) 
the imprinting systematically leads to the formation of 
a ^-soliton, i.e each shot produces a vortex ring and a 
vortex line. To lift this ambiguity, we establish, for each 
scenario, the constraints on the probabilities of (i) detect¬ 
ing no depletion (ii) detecting a horizontal vortex line (iii) 
detecting two nodal dots. Here, we focus on the four cen¬ 
tral slices near y = 0 (where horizontal vortex lines are 
observed). In the case of scenario (A): if p is the proba¬ 
bility to find a vortex line in one of the four slices, then 
the probability of detecting no depletion at all should be 
> 3p. Allowing for an inequality here accounts for the 
fact that one cannot exclude shots where the imprinting 
failed. This relation is true regardless of the probability 
of having a vortex ring. On the other hand, in the case 
of scenario (B): if we call p the probability of finding a 
vortex line in one of the four slices, then the probability 
of detecting two nodal points should be 3p. A statisti¬ 
cal analysis over ~ 90 slices recorded between t = 50 ms 
and t = 100 ms leads to the histogram shown in Fig. 
The probability of finding a line in one of the four cen¬ 
tral slices in this time interval is p = 27(4)% while the 
probabilities of detecting no depletion and a pair of nodal 
points are respectively p^ = I(I)% and pvR = 72(4)%. 
This result excludes scenario (A) and strongly supports 
scenario (B): the systematic formation of a 4>-soliton via 
a double puncture of the soliton plane. 

Generation of the histogram. The statistical sam¬ 
ple consists of 22 slices at p = —39 pm, 30 slices at 
y = —13 pm, 24 slices at p = 13 pm and 13 slices at 
p = 39 pm. The probability is obtained from several 
runs of random sampling of the data. In a given run, 
we randomly pick quadruplets of four slices (one slice 
at each position) until a set of 13 quadruplets is formed 
(13 is the smallest number of repetitions recorded for 
a given slice position). Out of a given set we extract 



no feature vortex line 

of points 


FIG. 5. Histogram of events between t=50ms and t=100ms. 
The probability of detecting no distinctive feature at all is 
P 0 = 1 ( 1 )%, the probability of detecting a horizontal line is p = 
27(4)% and the probability of detecting pairs of nodal points 
is pvR = 72(4)%. The shaded orange area is the interval 
consistent with 3p given the uncertainty on p. 


the probabilities of the three events: ’no feature’, 
’vortex line’ and ’pair of nodal points’. This proce¬ 
dure is repeated 6000 times, where each run provides 
a new set of 13 quadruplets and the corresponding 
probabilities p, pvR and p 0 . From the 6000 runs, we 
obtain the standard deviation of each probability which 
provides their interval of confidence (error bars in Fig.[^. 

Rapid-Ramp versus Time-of-Flight imaging. 

The images obtained after the rapid ramp show strong 
modulations of the density distribution. Based on these 
images, one might conclude that the phase imprinting 
leads to an intense disturbance in the density accompa¬ 
nied by shock waves. However, this is not the case. In 
order to clarify this point, we show, in Fig.[^ the time se¬ 
quence of the dynamics following the imprint by probing 
the density distribution after a time of flight expansion of 
4 ms, without any change of the scattering length. One 
observes only weak perturbations in the density profile. 
In particular, one can see two density wavefronts of small 
amplitude propagating in opposite direction towards the 
z = ±Rz edges of the cloud, that are identified as the fast 
wavefronts seen in the rapid ramp images. For a more 
direct comparison, the figure also shows the same time 
sequence recorded after time of flight combined with the 
slicing phase, which obviously suffers from a lower sig¬ 
nal to noise. This comparison provides evidence that the 
phase imprint performed in our experiment is indeed a 
gentle process, which predominantly affects the phase of 
the superfluid order parameter rather than the density. 
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FIG. 7. Snake instability of a planar dark soliton in a unitary Fermi gas, observed at early times after the phase imprint. The 
images represent a time series in time steps of 0.5 ms after the phase imprint. Top: optical density of the central slice near 
z = 0, after rapid ramp and time of flight. Bottom: detected minima (black) and reconstruction using Fourier analysis (gold). 
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FIG. 8. Tomographic images at various times t after the phase-imprint. The right column shows the structure of the solitary 
waves as if one would observe along the 2 ;-axis, reconstructed using the residuals of the tomographic images. 
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FIG. 9. Comparison of images with and without rapid ramp of the magnetic field to the BEG side of the Feshbach resonance. 
Top: images after 9 ms expansion with the rapid ramp as described previously. Bottom: images after 4 ms expansion without 
changing the Feshbach field. The expansion time is chosen to approximately match the radial size of the cloud with those 
obtained with rapid ramp. In each case we show (a) the integrated density (no slicing) and (b) the density of the central layer. 
Without rapid ramp, the amplitude of the initial wave fronts is strongly reduced. 
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Imprinting laser power 


FIG. 10. Integrated density after rapid ramp of a phase imprinted superfiuid with increasing imprinting power. As the power 
is increased (images from left to right) we observe the gradual change from a shallow and strongly bent soliton to a straight 
dark soliton of high contrast to two solitons. 







10 


Velocity Control of the Soliton. By increasing the 
intensity of the imprinting laser we observe (see Fig. 10) 
the gradual change from creating fast shallow solitons to 
slow high contrast solitons and, at even higher power, 
the formation of multiple solitons. The fast solitons 
are associated to a strong bending due to the radial 
inhomogeneity (see the discussion on the dominant 
snaking mode in the main text). As the imprinting 
laser intensity directly controls the phase jump across 
the soliton, this observation constitutes a qualitative 
demonstration of the soliton’s current-phase relation. 


Trapping potential. The trapping potential has the 
following form 


V{x,y,z) = + y^) 


+ Vo 


1 _ e w2 1 


+ may . 


( 1 ) 


The first term is a saddle point potential from the 
magnetic field curvature. The second term describes the 
optical dipole trap (ODT), with a trap depth V^jh — 
11 ( 1 ) kHz and a gaussian waist W = 125(5) yum. The 
ODT counters the anti-confinement in the radial x and 
y directions due to the magnetic saddle potential. The 
effective confinement in the x and y directions, in the ab¬ 
sence of gravity, can then be modeled as a harmonic po- 
tential with trapping frequency = ( ^^2 ~ — 

27r69(6)Hz. Since gravity is not completely compen¬ 
sated, the atoms experience an effective vertical accel¬ 
eration of a = 1.9(l)m/s^. This gradient also shifts the 
saddle point of the potential in the vertical direction by 
10 ym. 

Expanding the potential for small distances from the 
origin, as in Ref. m, one obtains 

V « ^m{cuy + + ujy) + Cy^ + Cx^y , ( 2 ) 

with ujx — — ^^±(1 — ^/3) and 

uOy — — 25 a;^(l — ^), where 5 = | ^^2 

and C = ^ are parameters quantifying the an- 

harmonicity of the potential. C is to be compared to 

muj^ 

Co = 2 ^’ where Ry = 61(9) ym is the Thomas-Fermi 
radius in the vertical direction. We note that the model 
of Ref. [19] does not include the correction to ujx and 
the Cx^y term. With the given definition for 6 and C 
and our experimental parameters, one has S — 0.019(4) 
and C = 0.16(8) Co. 


Dynamics during rapid ramp. The detection 
procedure requires a 9 ms interval between the time t 


of the dynamics that we want to observe and the time 
when the image is actually taken. During this phase, 
the cloud spends ~ 2 ms in the strongly interacting 



FIG. 11. a) Schematics of the rapid ramp. The solid line 
shows the ramp of the Feshbach field. At wait time t after 
the imprint, the field is ramped from 832G to 760G in 1 ms, 
and then ramped down to 576G over 1.5 ms. The field is 
held at that value for 4 ms, and then ramped back to 780G 
over 1 ms. The field is held at 780G for 1.5 ms before the 
spatially selective optical pumping light is applied and the 
cloud is imaged. The dashed line shows the corresponding 
ao/as, which is 0 at 832G, 4300“^ at 760G, 225“^ at 576G, 
and 6500“^ at 780G. b) Extraction of the effective time delay 
after the rapid ramp, from the extrapolation of the sound 
waves trajectories. Solid black circles: upper and lower sound 
wave positions as function of time t. The linear fit (solid gray 
lines) yields an effective lag of At 1.7 ms. 


regime (/ci?as)“^ ^ 1 where the expansion is hydrody¬ 
namic. For the rest of the rapid ramp, the cloud is so 
weakly interacting that the expansion can be essentially 
considered ballistic. Therefore, we estimate that the 
9 ms of the combined rapid ramp and expansion can 
be accounted for as an effective lag of ~ 2 ms in the 
dynamics of the observed images. One can notice this 
lag in Fig. showing a zoom on the early stage 

of the sound trajectories. By extrapolating the sound 
trajectory, we find an effective lag of 1.7ms, consistent 
with our estimate. 


Tracking the waves. In Fig. we show the dis¬ 
placement as a function of time for sound waves and de¬ 
fects, as extracted from the residuals in Fig. 4 of the main 
text (for t < 100 ms), as well as from the residual of in¬ 
tegrated images (for t > 100 ms). This figure provides an 
overview of the superfluid dynamics following the phase 
imprint. The various trajectories plotted here show the 
initial sound propagation, the second set of sound waves, 
the slowly moving soliton, progressively decaying into a 
single solitonic vortex, which then precesses in the su¬ 
perfluid. The precession is seen in the right panel as an 
oscillation along the ^-axis with a period of ~ 1.4 s. 
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FIG. 12. Displacement as a function of time for sound waves and defects. Left panel: location of the sound waves and the 
defect. This graph is obtained from the residuals displayed in Fig. 4 of the main text. Black solid circles: initial sound waves. 
Gray solid circles and band: second set of sound waves. Red circles: defect in the central axial cut. Blue squares: defect in the 
outer axial cuts. Right panel: precession motion (projected on the z-axis) of the solitonic vortex resulting from the cascade 
(red open circles). This graph is obtained from images of line of sight integrated densities. Red solid circles are the locations 
of the solitonic vortex as detected in the central axial cut of the central slice. 












